The role of chloroplast movement in C4 photosynthesis: a theoretical analysis using a three-dimensional reaction-diffusion model for maize.

Chloroplasts movement within mesophyll cells in C4 plants is hypothesized to enhance the CO2 concentrating mechanism, but this is difficult to verify experimentally. A three-dimensional (3D) leaf model can help analyse how chloroplast movement influences the operation of the CO2 concentrating mechanism. The first volumetric reaction-diffusion model of C4 photosynthesis that incorporates detailed 3D leaf anatomy, light propagation, ATP and NADPH production, and CO2, O2 and bicarbonate concentration driven by diffusional and assimilation/emission processes was developed. It was implemented for maize leaves to simulate various chloroplast movement scenarios within mesophyll cells: the movement of all mesophyll chloroplasts towards bundle sheath cells (aggregative movement) and movement of only those of interveinal mesophyll cells towards bundle sheath cells (avoidance movement). Light absorbed by bundle sheath chloroplasts relative to mesophyll chloroplasts increased in both cases. Avoidance movement decreased light absorption by mesophyll chloroplasts considerably. Consequently, total ATP and NADPH production and net photosynthetic rate increased for aggregative movement and decreased for avoidance movement compared with the default case of no chloroplast movement at high light intensities. Leakiness increased in both chloroplast movement scenarios due to the imbalance in energy production and demand in mesophyll and bundle sheath cells. These results suggest the need to design strategies for coordinated increases in electron transport and Rubisco activities for an efficient CO2 concentrating mechanism at very high light intensities.

In silico study of the role of cell growth factors in photosynthesis using a virtual leaf tissue generator coupled to a microscale photosynthesis gas exchange model.

Computational tools that allow in silico analysis of the role of cell growth and division on photosynthesis are scarce. We present a freely available tool that combines a virtual leaf tissue generator and a two-dimensional microscale model of gas transport during C3 photosynthesis. A total of 270 mesophyll geometries were generated with varying degrees of growth anisotropy, growth extent, and extent of schizogenous airspace formation in the palisade mesophyll. The anatomical properties of the virtual leaf tissue and microscopic cross-sections of actual leaf tissue of tomato (Solanum lycopersicum L.) were statistically compared. Model equations for transport of CO2 in the liquid phase of the leaf tissue were discretized over the geometries. The virtual leaf tissue generator produced a leaf anatomy of tomato that was statistically similar to real tomato leaf tissue. The response of photosynthesis to intercellular CO2 predicted by a model that used the virtual leaf tissue geometry compared well with measured values. The results indicate that the light-saturated rate of photosynthesis was influenced by interactive effects of extent and directionality of cell growth and degree of airspace formation through the exposed surface of mesophyll per leaf area. The tool could be used further in investigations of improving photosynthesis and gas exchange in relation to cell growth and leaf anatomy.

Using a reaction-diffusion model to estimate day respiration and reassimilation of (photo)respired CO2 in leaves.

Methods using gas exchange measurements to estimate respiration in the light (day respiration Rd ) make implicit assumptions about reassimilation of (photo)respired CO2 ; however, this reassimilation depends on the positions of mitochondria. We used a reaction-diffusion model without making these assumptions to analyse datasets on gas exchange, chlorophyll fluorescence and anatomy for tomato leaves. We investigated how Rd values obtained by the Kok and the Yin methods are affected by these assumptions and how those by the Laisk method are affected by the positions of mitochondria. The Kok method always underestimated Rd . Estimates of Rd by the Yin method and by the reaction-diffusion model agreed only for nonphotorespiratory conditions. Both the Yin and Kok methods ignore reassimilation of (photo)respired CO2 , and thus underestimated Rd for photorespiratory conditions, but this was less so in the Yin than in the Kok method. Estimates by the Laisk method were affected by assumed positions of mitochondria. It did not work if mitochondria were in the cytosol between the plasmamembrane and the chloroplast envelope. However, mitochondria were found to be most likely between the tonoplast and chloroplasts. Our reaction-diffusion model effectively estimates Rd , enlightens the dependence of Rd estimates on reassimilation and clarifies (dis)advantages of existing methods.

Localization of (photo)respiration and CO2 re-assimilation in tomato leaves investigated with a reaction-diffusion model.

The rate of photosynthesis depends on the CO2 partial pressure near Rubisco, Cc, which is commonly calculated by models using the overall mesophyll resistance. Such models do not explain the difference between the CO2 level in the intercellular air space and Cc mechanistically. This problem can be overcome by reaction-diffusion models for CO2 transport, production and fixation in leaves. However, most reaction-diffusion models are complex and unattractive for procedures that require a large number of runs, like parameter optimisation. This study provides a simpler reaction-diffusion model. It is parameterized by both leaf physiological and leaf anatomical data. The anatomical data consisted of the thickness of the cell wall, cytosol and stroma, and the area ratios of mesophyll exposed to the intercellular air space to leaf surfaces and exposed chloroplast to exposed mesophyll surfaces. The model was used directly to estimate photosynthetic parameters from a subset of the measured light and CO2 response curves; the remaining data were used for validation. The model predicted light and CO2 response curves reasonably well for 15 days old tomato (cv. Admiro) leaves, if (photo)respiratory CO2 release was assumed to take place in the inner cytosol or in the gaps between the chloroplasts. The model was also used to calculate the fraction of CO2 produced by (photo)respiration that is re-assimilated in the stroma, and this fraction ranged from 56 to 76%. In future research, the model should be further validated to better understand how the re-assimilation of (photo)respired CO2 is affected by environmental conditions and physiological parameters.

Impact of anatomical traits of maize (Zea mays L.) leaf as affected by nitrogen supply and leaf age on bundle sheath conductance.

The mechanism of photosynthesis in C4 crops depends on the archetypal Kranz-anatomy. To examine how the leaf anatomy, as altered by nitrogen supply and leaf age, affects the bundle sheath conductance (gbs), maize (Zea mays L.) plants were grown under three contrasting nitrogen levels. Combined gas exchange and chlorophyll fluorescence measurements were done on fully grown leaves at two leaf ages. The measured data were analysed using a biochemical model of C4 photosynthesis to estimate gbs. The leaf microstructure and ultrastructure were quantified using images obtained from micro-computed tomography and microscopy. There was a strong positive correlation between gbs and leaf nitrogen content (LNC) while old leaves had lower gbs than young leaves. Leaf thickness, bundle sheath cell wall thickness and surface area of bundle sheath cells per unit leaf area (Sb) correlated well with gbs although they were not significantly affected by LNC. As a result, the increase of gbs with LNC was little explained by the alteration of leaf anatomy. In contrast, the combined effect of LNC and leaf age on Sb was responsible for differences in gbs between young leaves and old leaves. Future investigations should consider changes at the level of plasmodesmata and membranes along the CO2 leakage pathway to unravel LNC and age effects further.

Mesophyll conductance and reaction-diffusion models for CO2 transport in C3 leaves; needs, opportunities and challenges.

One way to increase potential crop yield could be increasing mesophyll conductance gm. This variable determines the difference between the CO2 partial pressure in the intercellular air spaces (Ci) and that near Rubisco (Cc). Various methods can determine gm from gas exchange measurements, often combined with measurements of chlorophyll fluorescence or carbon isotope discrimination. gm lumps all biochemical and physical factors that cause the difference between Cc and Ci. gm appears to vary with Ci. This variability indicates that gm does not satisfy the physical definition of a conductance according to Fick's first law and is thus an apparent parameter. Uncertainty about the mechanisms that determine gm can be limited to some extent by using analytical models that partition gm into separate conductances. Such models are still only capable of describing the CO2 diffusion pathway to a limited extent, as they make implicit assumptions about the position of mitochondria in the cells, which affect the re-assimilation of (photo)respired CO2. Alternatively, reaction-diffusion models may be used. Rather than quantifying gm, these models explicitly account for factors that affect the efficiency of CO2 transport in the mesophyll. These models provide a better mechanistic description of the CO2 diffusion pathways than mesophyll conductance models. Therefore, we argue that reaction-diffusion models should be used as an alternative to mesophyll conductance models, in case the aim of such a study is to identify traits that can be improved to increase gm.

A two-dimensional microscale model of gas exchange during photosynthesis in maize (Zea mays L.) leaves.

CO2 exchange in leaves of maize (Zea mays L.) was examined using a microscale model of combined gas diffusion and C4 photosynthesis kinetics at the leaf tissue level. Based on a generalized scheme of photosynthesis in NADP-malic enzyme type C4 plants, the model accounted for CO2 diffusion in a leaf tissue, CO2 hydration and assimilation in mesophyll cells, CO2 release from decarboxylation of C4 acids, CO2 fixation in bundle sheath cells and CO2 retro-diffusion from bundle sheath cells. The transport equations were solved over a realistic 2-D geometry of the Kranz anatomy obtained from light microscopy images. The predicted responses of photosynthesis rate to changes in ambient CO2 and irradiance compared well with those obtained from gas exchange measurements. A sensitivity analysis showed that the CO2 permeability of the mesophyll-bundle sheath and airspace-mesophyll interfaces strongly affected the rate of photosynthesis and bundle sheath conductance. Carbonic anhydrase influenced the rate of photosynthesis, especially at low intercellular CO2 levels. In addition, the suberin layer at the exposed surface of the bundle sheath cells was found beneficial in reducing the retro-diffusion. The model may serve as a tool to investigate CO2 diffusion further in relation to the Kranz anatomy in C4 plants.

Three-dimensional microscale modelling of CO2 transport and light propagation in tomato leaves enlightens photosynthesis.

We present a combined three-dimensional (3-D) model of light propagation, CO2 diffusion and photosynthesis in tomato (Solanum lycopersicum L.) leaves. The model incorporates a geometrical representation of the actual leaf microstructure that we obtained with synchrotron radiation X-ray laminography, and was evaluated using measurements of gas exchange and leaf optical properties. The combination of the 3-D microstructure of leaf tissue and chloroplast movement induced by changes in light intensity affects the simulated CO2 transport within the leaf. The model predicts extensive reassimilation of CO2 produced by respiration and photorespiration. Simulations also suggest that carbonic anhydrase could enhance photosynthesis at low CO2 levels but had little impact on photosynthesis at high CO2 levels. The model confirms that scaling of photosynthetic capacity with absorbed light would improve efficiency of CO2 fixation in the leaf, especially at low light intensity.

Modelling the relationship between CO2 assimilation and leaf anatomical properties in tomato leaves.

The CO2 concentration near Rubisco and, therefore, the rate of CO2 assimilation, is influenced by both leaf anatomical factors and biochemical processes. Leaf anatomical structures act as physical barriers for CO2 transport. Biochemical processes add or remove CO2 along its diffusion pathway through mesophyll. We combined a model that quantifies the diffusive resistance for CO2 using anatomical properties, a model that partitions this resistance and an extended version of the Farquhar-von Caemmerer-Berry model. We parametrized the model by gas exchange, chlorophyll fluorescence and leaf anatomical measurements from three tomato cultivars. There was generally a good agreement between the predicted and measured light and CO2 response curves. We did a sensitivity analysis to assess how the rate of CO2 assimilation responds to changes in various leaf anatomical properties. Next, we conducted a similar analysis for assumed diffusive properties and curvature factors. Some variables (diffusion pathway length in stroma, diffusion coefficient of the stroma, curvature factors) substantially affected the predicted CO2 assimilation. We recommend more research on the measurements of these variables and on the development of 2-D and 3-D gas diffusion models, since these do not require the diffusion pathway length in the stroma as predefined parameter.